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摘 要 : 针对 全 局 优化 问题 , 基于 一 类 支持 向 量 数据 描述 (SVDD) 和 已 有 的 根系 生长 算法 提出 一 种 新 的 智能 优化 算法 
一 一 根 生 群 优化 算法 ,将 根系 划分 为 主根 群体 和 侧根 群体 。 基 于 SVDD 描述 主根 群体 的 生长 行为 ， 将 土壤 中 养分 浓度 
最 高 的 位 置 作为 全 局 优化 的 目标 ， 构 建 了 根系 生长 模型 ; 分 析 了 RGSO 的 数学 模型 ， 从 理论 上 证 明了 RGSO 的 收敛 
性 。 在 实验 中 ， 与 当前 最 先进 的 其 他 三 个 算法 进行 综合 比较 ， 并 观察 了 不 同 参数 对 优化 效果 的 影响 。 实 验 结 果 验 证 了 
RGSO 的 收敛 性 和 有 效 性 ， 表 明 RGSO 是 一 种 解决 全 局 优化 问题 的 有 效 算 法 。 
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Abstract: Aiming at the global optimization problem, this paper proposed a new intelligent algorithm called RGSO, based on 
SVDD and the root algorithm existed. The root system differentiated into the taproot and the lateral root group. Inspired by this 
growth behavior, this paper described the growth behavior of the taproot tips based on SVDD. The algorithm constructed the 
root growth model, and used the point of the highest concentration of nutrients in the soil as the target of global optimization. 
This paper analyzed the mathematical model of RGSO, and proved its convergence theoretically. In the experiment, this paper 
compared RGSO with the other three advanced algorithms, and tested the optimization effect with different parameters. The 


result of the experiment verified the convergence and effectiveness of RGSO, and it indicates that RGSO is an effective algorithm 


to solve global optimization problem. 


— Key Words: root algorithm; clustering mechanism; SVDD; optimization problem 


其 中 : x An 维 向 量 , /为 目标 函数 ，h 和 g 为 约束 条 件 。 
生物 智能 算法 是 一 种 新 兴 的 演化 计算 技术 ， 其 研究 对 象 包 
自然 生物 现象 遵循 着 某 些 内 在 规律 止 。 几 个 世纪 以 来 ， 科 括 生 物 自身 的 行为 名 以 及 生物 所 处 的 自然 环境 ， 目 的 是 找 出 其 
学 家 们 试图 识别 和 记录 自然 现象 背后 的 物理 定律 中。 受 各 种 自 中 的 数学 规律 。 智 能 优化 算法 要 解决 的 一 般 是 最 优化 问题 。 传 
然 生物 现 象 的 启发 ， 计 算 机 工作 者 们 以 自然 界 生物 的 行为 B 为 。” 统 的 智能 算法 都 具有 自 适 应 的 特点 ， 如 遗传 算法 四、 免疫 算法 


0 引言 


Ds 


基础 ， 抽 取出 用 于 求解 大 规模 复杂 问题 的 计算 技术 和 方法 。 最 。 中 、 粒 子 群 算法 外 等 。 
无 化 问题 外 是 一 个 极其 古老 的 课题 。 目 前 ， 优 化 问题 在 工程 技 粒子 群 优 化 算法 (PSO) 是 一 种 源 于 对 鸟 群 捕食 行为 的 研究 
术 、 经 济 管理 、 交 通 运输 、 系 统 控制 、 人 工 智 能 、 模 式 识别 等 ” 而 发 明 的 进化 计算 技术 ,最 先 由 Barnhart 博士 和 Kennedy 博士 
领域 得 到 了 成 功 的 推广 和 应 用 。 提出 四。2014 年 ，Zhang 等 人 09 提 出 了 原始 的 根 生 算法 ， 模 拟 
最 优化 问题 的 数学 描述 如 下 : 了 植物 根系 的 生长 ， 将 寻找 土壤 中 养分 浓度 最 高 的 点 作为 寻 优 
h(x) =0; 的 目标 。 支 持 向 量 数据 描述 (SVDD) 是 用 于 一 类 分 类 问题 的 分 类 
求解 min f(x), s.t. 
g(x) 20 技术 。2013 4F, Niazmardi 等 人 0 提出 了 一 个 改进 的 聚 类 框架 
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作为 在 聚 类 时 缺乏 足够 高 的 高 品质 的 训练 数据 的 解决 方案 。 

些 基 于 传统 算法 的 改进 算法 也 被 不 断 提 出 。 文 献 [12] 基 
于 前 馈 神 经 网 络 的 反 向 传播 学 习 ， 用 于 优化 FPGA 芯片 。 文 献 
[13] 提 出 带 干 扰 抵消 的 协同 传输 的 蚁 群 算法 。 文 献 [14] 基 于 多 
标的 模拟 退火 和 遗传 算法 ， 用 于 优化 三 相 高 温 超 导 变 压 器 。 文 
献 [15] 使 上 
文献 [16] 提 出 了 一 个 能 够 在 动态 环境 中 解决 时 变 非 线性 问题 的 
人 工 免 疫 算法 。 文 献 [17] 通 过 强化 单个 粒子 的 搜索 能 力 提 出 了 


一 个 自主 学 习 粒 子 群 优化 算法 。 文 献 [18] 对 粒子 群 算法 的 邻居 
结构 个 具有 综合 学 习 能 力 的 粒子 群 


进行 了 改进 ， 从 而 提出 了 
多 化 算法 。 近 年 来 ,新 提出 的 智能 算法 不 断 地 涌现 。 文 献 [19] 基 
于 改进 的 社会 蜘蛛 优化 算法 应 用 于 电磁 优化 ， 文 献 [20] 提 出 绒 


包 菌 算法 用 于 解决 斯 坦 纳 树 问题 ， 文 献 [21] 提 出 蚊子 追踪 算法 


基于 以 上 相关 工作 ， 本 文 所 设计 根 生 优化 模型 ， 规 定 了 算 
法 的 自 适 应 学 习 方 式 。 通 过 自 适应 学 习 的 特性 ， 算 法 达到 了 全 
局 优化 的 目的 。 大 多 数 自然 启发 算法 都 有 着 一 定 程度 的 有 效 性 ， 
但 是 也 有 着 相同 的 缺点 ， 即 收敛 速度 较 慢 ， 容 易 陷入 早熟 和 局 
部 最 优 。 本 文 所 提出 的 基于 SVDD 的 根 生 群 优化 算法 ， 即 是 基 
于 己 提 出 的 根 生 算 法 ， 针 对 上 述 的 缺点 ， 对 原 算法 的 改进 。 
RGSO 刻画 了 根系 群体 的 生长 方式 以 及 个 体 间 的 相互 作用 ， 与 
其 他 算法 相 比 , 具有 更 高 的 自 适 应 和 自我 学 习 能 力 . 因 此 RGSO 
在 求解 大 规模 复杂 问题 时 ， 具 有 更 高 的 有 效 性 。 


1 ”根系 群体 和 生长 行为 
根系 群体 


1.1 


精英 非 支 配 排 序 遗 传 算法 解决 多 目标 路 径 规 划 问 题 。 


向 着 这 个 更 好 的 位 置 生长 。 


[ 
适应 度 


DA: 


值 ， 个 


主根 根 尖 生 长 策 


体 适 应 度 高 于 该 平 


均 


于 种 群 的 多 样 性 ， 每 个 根 尖 所 处 的 位 置 不 同 ， 其 相应 的 


也 不 相同 。 依 据 种 群 中 根 尖 的 适应 度 ， 
名 和 全 
根 的 依据 是 个 体 适 应 度 的 大 小 ， 计 算 晶 


WAR ARR AE ARE 


值 的 根 尖 ， 将 被 


根 尖 的 生长 策略 
各 。 划 分 主根 与 侧 


所 有 个 体 的 适应 度 平均 


作为 主根 根 尖 ， 个 


体 适应 度 低 于 该 平均 值 的 生长 点 ， 将 被 作为 侧根 生长 点 。 


图 1 植物 根系 生长 形态 


1.4 生长 素 因子 
物 根系 在 正常 的 生长 状态 下 ， 会 产生 一 种 生长 激素 以 便 


植 
促进 根 


的 适应 值 水 平 。 当 前 根 尖 所 处 位 置 


优秀 ( 


1.5 土壤 养分 区 域 分 布 模型 


域 集中 
那么 即 
了 确定 


尖 的 生长 。 生 长 素 


越 小 ) ， 则 生长 素 


际 上 ， 土 壤 中 养分 
的 分 布 。 一 般 的 ,如 


大 


的 


H, 
N 


的 大 小 反映 了 当前 根 尖 在 当前 群体 中 


的 适应 值 相 较 于 其 他 个 体 越 


子 相 较 于 其 他 个 体 越 大 。 


分 布 不 是 均匀 的 ， 总 体 上 呈现 出 区 


土壤 中 多 个 位 


的 养分 浓度 较 高 ， 


认为 


植物 根系 是 由 一 系列 的 根 尖 构成 ， 根 尖 是 根 中 生命 活动 最 
活跃 的 部 分 ， 也 是 根 的 生长 、 延 长 及 吸收 水 分 的 主要 部 分 ， 因 


建立 土 


这 多 个 位 置 所 确定 的 空间 
这 一 区 域 的 位 置 , 引入 一 类 支持 


廊 向 量 空间 模型 ， 即 建立 一 个 


将 待定 


此 根 的 生长 依赖 于 根 尖 部 位 的 生长 。 根 尖 的 生长 过 程 与 其 生长 


环境 息息相关 ， 最 主要 的 影响 因素 即 是 土壤 中 养分 的 浓度 。 事 
实证 明 , 植 物 的 根 尖 总 是 向 着 土壤 中 养分 浓度 最 高 的 地 方 生长 ， 
如 图 1 所 示 。 
1.2 根系 初始 状态 

植物 的 根系 是 由 一 粒 种 子 发 育 而 来 ， 即 根系 群体 在 初始 状 
态 下 根 尖 的 个 体 数 为 1。 


1.3 主根 子 群 和 侧根 子 群 
在 RGSO 中 ,， 根 生 群 体 被 划分 为 两 个 子 群 体 ， 即 主根 子 群 
和 侧根 子 群 。 以 下 是 两 个 子 群 的 定义 。 
定义 1 主根 子 群 。 胚 根 细胞 分 裂 和 伸 长 向 下 垂直 生长 的 
根 , 是 植物 体 最 早出 现 的 根 ,也 称 直 根 或 初生 根 。 大 多 数 裸子 植 
物 和 双子 叶 植物 主根 继续 生长 ， 明 显 而 发 达 。 植 物 主要 依靠 主 
根 的 不 断 向 下 生长 ， 寻 找 土壤 中 养分 浓度 最 高 的 位 置 。 
定义 2 ”侧根 子 群 。 土 壤 中 的 主根 不 断 地 生长 ， 长 到 一 定 
的 程度 便 伸 出 侧 支 ， 形 成 庞大 的 根系 。 形 成 侧 支 的 主要 动因 是 
为 了 更 大 范围 地 搜索 土壤 中 浓度 较 高 的 位 置 ， 一 旦 发 现 养分 浓 
度 更 高 的 位 置 ， 就 会 将 这 一 情况 反馈 给 根系 ， 整 个 根系 也 随 之 


这 一 区 
1.6 

根 
生长 行 


pa 
EX 


的 感知 。 
段 设 2 ” 根 尖 在 和 9 


Fy 
{EL 


p: 
EX 


或 。 


生长 行为 
系 的 生长 行为 主要 包括 根 尖 的 分 支 和 生长 ， 现 将 根系 的 


区 域 的 多 个 向 量 包 库 ， 得 到 其 


为 作 以 下 的 限定 : 


设 1 根 尖 的 生长 方向 主要 依靠 


区 域 养分 浓度 较 高 。 为 
向 量 数据 描述 (SVDD)， 
尽 可 能 小 的 向 量 空间 球 ， 
球 心 位 置 ， 以 该 球 心 表征 


ai 


则 主根 会 进行 分 支 。 


体能 够 指引 落后 个 体 的 生长 。 


较 差 的 


根 尖 。 


1.6.1 主根 的 分 支 


根 
群体 中 


Zt 
根 尖 数目 


Ri 4 ” 当 根 尖 总 数 达到 一 定 的 阐 值 后 ， 淘 汰 一 些 适 应 


自身 对 土壤 中 养分 浓 


Xe 


长 的 过 程 中 ， 如 果 群 体 根 尖 总 数 未 达到 


设 3 不同 的 根 尖 个 体 之 间 的 行为 会 相互 影响 ， 优 秀 个 


KE 


一 粒 种 子 〈 一 个 根 尖 ) 发 育 而 来 ， 根 系 的 生长 即 为 
的 不 断 扩充 。 在 根系 的 生长 过 程 中 ， 优 秀 的 个 


体会 在 其 生长 空间 周围 进行 分 支 ， 根 尖 的 分 支 会 造成 群体 根 尖 


总 数 的 
且 每 次 


增加 。 在 本 文中 , 限定 只 有 主根 根 尖 个 体 才 能 进行 分 支 ， 


分 支 的 数目 


恒定 。 
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1.6.2 主根 的 生长 
主根 根 尖 的 生长 较 慢 ， 受 两 个 因素 影响 : 
a) 角 度 因素 。 由 于 主根 总 是 向 着 土壤 中 养分 浓度 较 高 的 位 
置 生长 ,因此 ,角度 的 朝向 是 当前 养分 浓度 较 高 的 区 域 的 位 置 ， 
该 区 域 是 由 所 有 优 于 当前 个 体 的 根 尖 所 确定 的 空间 构成 ， 该 空 
可 以 利用 SVDD 确定 ， 并 且 求 出 其 球 心 位 置 ， 则 当前 根 尖 应 
该 朝 着 该 球 心 生 长 。 
b) 长 度 因素 。 该 长 度 由 种 群 所 有 个 体位 置 的 ] 
个 体位 置 的 距离 和 当前 个 体 的 生长 素 因 子 决定 。 
1.6.3 侧根 的 生长 
侧根 根 尖 的 生长 较 快 ,一 般 认为 ,侧根 根 尖 的 生长 较 散乱 ， 
生长 角度 与 方向 都 较为 随机 。 
1.7 淘汰 机 制 
事实 上 ， 根 系 群体 的 根 尖 数目 是 有 上 限 的 ， 这 个 上 限 ， 也 
被 称 作 种 群 个 体 数目 的 平衡 点 。 当 总 群 个 体 数目 超过 这 个 平衡 
点 后 ， 种 群 会 淘汰 那些 落后 的 个 体 ， 以 达到 整体 个 数 的 平衡 。 
淘汰 的 准则 是 剔除 那些 适应 值 最 差 的 个 体 ， 直 到 整体 个 数 重 新 
达到 平衡 点 。 


2 ， 根 生 群 优化 数学 模型 和 算法 
根据 根系 的 形成 过 程 、 根 系 的 分 群 过 程 以 及 根 尖 的 分 支 和 


生长 行为 ， 得 到 流程 图 (图 2) 。 


初始 化 种 子 的 位 置 以 及 各 参数 * 


计算 个 体 的 适应 度 及 生长 素 浓度 。 


划分 主根 子 群 和 侧根 子 群 
及 主根 根 尖 进行 分 支 。 


可 


匀 值 点 到 当前 


= 


主根 和 侧根 依据 不 同 的 策略 进行 生长 


个 体 数 达到 。 
RAY 


>| mam ERENT. 


是 否 达到 售 
止 准则 J 


图 2 根 生 群 优化 生长 流程 


2.1 根系 群体 的 数学 模型 
定义 3 根系 群体 的 初始 状态 是 一 颗 种 子 ， 即 算法 在 初始 
状态 ， 只 有 一 个 根 尖 。 这 个 根 尖 的 位 置 是 在 可 行 域 中 随机 生成 
的 ， 整 个 根系 的 生成 都 依赖 于 这 个 根 尖 。 

who=l+A:(u—)) (1) 


其 中 : 和 /分 别 表示 n 维 的 可 行 域 上 下 界 向 量 ,，4 是 一 个 [0, 1] 


的 随机 变量 。 


根据 根系 子 群 的 划分 ， 根 系 群体 被 分 为 主根 子 群 和 侧根 子 


群 。 


定义 4 对 于 任意 的 根 尖 个 体 Ri:，Rie S,，S = {M,L}, 其 
H, S 为 根 尖 群体 集合 ，M 为 主根 子 群集 合 ， 为 侧根 子 群集 


A 
FI o 


算法 1 根系 子 群 分 类 


HA: 适应 度 函数 fit(x) 以 及 各 参数 


1: 计 算 所 有 根 尖 的 适应 度 fit(xi) 
2: 计 算 适应 度 平均 值 aver = > fit) 
3: for each root tips 
4: if fit(xi) 2aver 
5: ”该 根 尖 属于 主根 
6: else 
7: 该 根 尖 属 于 侧根 
8: end if 
9: end for 
输出 ， 根 尖 个 体 分 类 
2.2 群生 长 算法 
2.2.1 新 的 根 尖 位 置 


群生 长 包括 新 的 根 尖 位 置 的 


长 。 
定义 5 ”新 的 根 尖 的 形成 位 置 与 原 根 尖 位 
点 的 形成 位 置 定义 如 下 : 
Xew ZX, 


nev 


成 ， 根 尖 分 支 以 及 


相关 


„a + AX 


其 中 :xow 表 示 旧 的 根 尖 的 位 置 , Ax 是 


布 变量 。 

2.2.2 主根 分 支 
定义 6 REE 
算法 2. 根 尖 分 支 算法 


输入 : 根 尖 个 体 的 位 置 


1: for each root tip xi 


2: if xiem 


3: 根据 式 (2) 产 生 新 的 根 尖 位 置 


4: end if 


3: end for 


输出 ;新 产生 的 根 尖 位 置 


个 服从 N(0,1) 


根 才 进行 分 支 ， 侧 根 不 进行 分 支 。 


2.2.3 基于 SVDD 的 土壤 区 域 分 布 


定义 7 优 于 当前 根 尖 的 所 


Be 


了 根 尖 个 体 所 确定 的 


(ROE) =R +C 2 


(x -0V (x 


OVER +E 


根 尖 的 生 


， 新 生长 


(2) 
的 正 态 分 


区 域 由 包 


这 些 个 体 向 量 的 空间 超 球 确定 ， 该 超 球 的 球 心 设 为 0， 半 径 
R 


(3) 


其 中 :xie 4, 4 是 所 有 优 于 当前 个 体 的 根 尖 集合 , O 是 超 球 的 球 


心 ,R 是 超 球 的 半径 ,是 松弛 变 直 


a,C 是 惩罚 因子 ,基于 


F SVDD, 
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可 以 根据 4 计算 出 O 和 RR。 
2.2.4 主根 和 侧根 生长 策略 
在 RGSO 中 , 根系 群体 分 为 主根 子 群 和 侧根 子 群 
群 分 别 拥 有 不 同 的 生长 方式 。 

对 于 主根 子 群 ， 其 生长 较为 缓慢 ， 每 次 生长 都 要 将 根系 群 
体 中 较为 优秀 的 个 体 的 信息 作为 参考 ， 调 整 自身 的 生长 幅度 。 
较为 优秀 的 个 体 ， 即 为 适应 值 小 于 当前 个 体 的 所 有 其 他 个 体 。 


两 种 子 


= 


对 于 侧根 子 群 ， 其 生长 较为 迅速 ， 侧 根 的 生长 依照 随机 游 走 的 
原则 。 
定义 8 主根 个 体 每 次 生长 的 角度 由 其 他 优 于 当前 个 体 的 
根 尖 决 定 。 
O, — x; 
“jo -x] (4) 


其 中 :Oi 是 所 有 优 于 xi 的 个 体 所 确定 的 超 球 球 心 。 
定义 9 生长 素 反 映 了 根 尖 个 体 当 前 的 适应 值 水 平 ， 适 应 
值 越 小 ， 生 长 素 值 越 大 。 


fx) ) 
DSa) 6) 


jes 


E=exp(- 


其 中 : S 表示 当前 群体 的 所 有 根 尖 个 体 集合 。 
定义 10 主根 根 尖 每 次 生长 长 度 y 与 当前 的 生长 素 值 相 
关 ， 生 长 素 值 越 大 ， 则 需要 生长 的 长 度 越 小 ， 反 之 ， 则 越 大 。 


y= Brand E r -a (6) 

其 中 ,nn 是 当前 个 体 总 数 ，E 是 生长 素 因子 ，rand 是 一 个 [0, 1] 
的 随机 数 。 

定义 11 缩放 算 子 o 与 当前 根 尖 的 适应 值 相关 ， 适 应 值 


越 小 ， 说 明 当 前 根 尖 越 优 秀 ， 则 保留 本 来 位 置 的 可 能 越 大 。 


n- L10 


1 (7) 
2 fa 
定义 12 主根 根 尖 的 位 置 更 新 ， 定 义 如 下 : 
X(t+N=a-x,(tht+a-0-v (8) 


RE, ORRE i 个 主根 根 尖 在 1 时 刻 的 位 置 ,a 是 控制 参数 ， 
取 值 为 [0, 1]，0 表示 生长 角度 ，y 表示 生长 的 长 度 。 
定义 13 ”侧根 根 尖 的 位 置 更 新 
x(t +1) = rand - B-x,(t) (9) 
RE, ORRE i 个 侧根 根 尖 在 1 时 刻 的 位 置 , 采用 随机 游 走 
的 生长 策略 ，p 是 控制 参数 ，rand 是 一 个 [0, 1] 的 随机 数 。 
算法 3. 根 尖 生长 算法 


输入 : 根 尖 个 体 的 位 置 、 类 别 以 及 各 参数 


1: for each root tip 

2: if xi eM 

3: 根据 式 (8) 计 算 主根 位 置 
4: else 


5: 根据 式 (9) 计 算 侧根 位 置 


6: end if 
7: end for 
输出 ， 更 新 后 的 根 尖 位 置 
种 群 的 每 一 次 分 支 的 过 程 都 会 使 得 种 群 根 尖 数 增加 ， 
本 文 设 


KIE, 


了 种 群 个 体 总 数 的 上 限 Num。 当 种 群生 长 点 总 数 超过 


Num 时 ， 则 将 会 依据 每 个 生长 点 的 适应 度 对 种 群 中 处 于 弱势 的 
个 体 进行 淘汰 。 淘 汰 的 准则 是 ， 每 次 淘汰 种 群 中 适应 度 处 于 最 


AMMA, BBS BUN F Bere BE Num. 
3 ” 根 生 优化 模型 的 理论 分 析 


3.1 ”收敛 性 证 明 
引 理 1 在 等 比 数列 {a} 
收敛 的 充分 条 件 是 |g |<1。 


引 理 2 一 阶 常 系数 线性 非 齐 次 方程 : 
x! + ax* = f (x) 

解 的 收敛 条 件 是 加 <1 

证 明 


上 述 差 分 方程 对 应 的 特征 方程 是 1~-a=0， 解 得 14= 
则 该 方程 对 应 的 齐 次 方程 通 解 为 x(k) = Cat 


中 ， 通 项 an ag", a0, MW fan} 


a 


为 保证 方程 的 解 收敛 ， 必 须 保证 齐 次 通 解 的 收敛 ， 即 


， k 
lima =0， 有 lol<1 


证 毕 。 


定理 1 对 于 向 量 集合 4， 其 中 每 一 个 向 量 都 是 n 维 ， 必 
然 能 够 找到 一 个 n 维 的 超 球面 ， 将 集合 A PETA Te 
定理 2 车 0<a<1,0<b<1， 则 有 
0<ļa-b|<1 
定理 3 根 生 群 优化 算法 的 主根 生长 行为 是 收敛 的 。 
证 明 
为 简化 证 明 过 程 ， 对 于 本 文中 的 n 维 向 量 ， 只 考察 其 在 1 


AS 


EE 时 的 情况 ， 即 n=1。 
主根 的 生长 表达 式 为 


x(t+D)=0-x(t)+a-0-y 


将 0 展开 得 
ay 
|O,(t) - x,(1)| 


将 上 式 转换 为 差分 方程 的 形式 : 


x(t+l =a-x,()4 (O,(t) - x,(t)) 


a-y æ- y- Olt) 
t +1) +(————— - 0): x(t) = 
neva EO ae |O(t) -x)| 
ay 
Sid 上 = 一 一 一 一 一 
BAS om-x 
重 写 差 分 方程 得 
a:-y:-O(t) 
t+1)+k-x(t) = 一 -一 一 一 一 
x(t +1) +k- x(t) ow) =x) 


He kr y 展开 得 
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k= 


a-(1—E)-rand - 


L Fs0-x0 


(a) 
|O(t) - x(t) 


1 
n LO — x(t) 


& p= 


O(t) - x(t) 


WW k =a- (1 — E) - rand - p —- œ 
4 t> Dy, A 


lim p =lim 


too tra 


1 
pe — x(t) 
O(t) - x(t) 


MO<a<1, 0<1-E<1, 0<rand<1, WA 
0<a-(1—-E)-rand:p<1, 


又 0<o<1, 根 据 定 理 2, 1174 0<|a-(1-E)-rand- p-oa|<1, 


即 O<|k\<1, 根据 引 理 2 可 知 ， 差 分 方程 的 解 收敛 的 充 要 条 件 即 


是 0<|k| <1。 


证 毕 。 


3.2 


在 实验 中 可 以 证 


有 效 性 分 析 
于 植物 根系 的 向 营养 性 ， 植 物 的 根 尖 总 是 朝 着 土壤 中 养 
分 浓度 较 高 的 区 域 生 长 ，RGSO 很 好 地 模拟 了 这 一 特性 。 
体现 种 群 的 多 样 性 ，RGSO 将 生长 点 分 为 主根 生长 点 和 侧根 生 
长 点 ,根据 差分 方程 的 收敛 性 , 可 以 推导 出 寻 优 结果 的 收敛 性 。 


H 


为 了 


H, RGSO 确实 可 以 快速 有 效 地 找到 极 值 点 。 


4 ”实验 及 结果 分 析 


实验 环境 : 
Shenteng 6800， 该 集群 拥有 8 个 计算 节点 和 一 个 总 控 节 点 。 每 


一 个 计算 节点 是 一 台 高 性 能 服务 器 ， 内 存 为 24 GB, 
GHz 中 央 处 理 器 。 所 


Linux 7. 
4.1 基准 函数 
实验 中 | 


于 测试 的 


MATLAB R2014b on the HPCC of Lenovo 


四 核 2.4 
服务 器 的 操作 系统 是 Red Hat Enterprise 


车 准 函数 包括 10 个 在 智能 进化 算法 中 


经 常用 到 的 基 7 


住 函数 ， 这 些 基准 函数 如 表 2 所 示 。 
42 生长 策略 与 参数 分 析 


实验 中 的 迭代 总 数 设置 为 1000 次 ， 根 生 优 化 算法 的 各 参 


数 设 


如 表 1 所 示 ， 


个 参数 的 设 定 都 是 经 过 多 次 实验 的 测 


试 以 确保 算法 的 效果 达到 最 好 。 


主根 的 生长 是 将 群体 中 较为 优秀 的 个 
一 个 包含 当前 所 有 较 优秀 个 体 的 超 


体 作 为 参考 ， 即 建立 


球 ， 以 该 超 球 球 心 作为 当前 


个 体 追 随 的 目标 ， 
程度 跳出 局 部 最 优 的 作用 。 侧 根 生长 依照 随机 游 走 的 原则 ， 即 


使 在 主根 生长 陷入 局 部 最 优 时 ， 
参数 a 起 到 调节 主根 生长 步 长 的 作用 ， 


目的 是 最 大 化 利 息 ， 起 到 一 定 


所 有 有 用 


期 望 能 够 跳出 局 部 最 优 。 


AE HAT 
MARAI 法 


算法 的 数值 寻 优 结果 随 a 值得 变化 如 图 
数值 取 以 


3 Arm, B 
自然 常数 e 为 底 的 对 数 。a 在 [0.5, 


中 纵 轴 的 
取 值 时 ， 算 


0.6] 


过 小 ， 则 算法 精度 较 差 ，b 
起 到 调节 侧根 生长 步 长 的 作 
时 ， 配 合 以 上 参数 ， 算 法 的 效果 达到 最 好 。 


个 体 适应 度 (Sphere) 


法 性 能 达到 最 好 。 参 数 b 影响 算法 的 收敛 速度 和 精度 , 若 b 值 
ITA, WEE RE. BR 8 
1, ABM, p 取 值 为 [5, 10] 


Fitness Value(log) 


Iteration Number 


#1 RGSO 参数 设 定 


1 1 h 1 1 L 1 1 1 J 
0 100 200 300 400 500 600 700 800 900 1000 


图 3 在 不 同 w 值 下 算法 寻 优 结果 (Sphere, D=60) 


参数 fe 意义 

r 1 根 尖 初 始 数量 
Num 100 群体 规模 

b 3 分 支 数 

a [0, 1] 参数 

B [5,10] 参数 


4.3 与 其 他 算法 比较 及 分 析 


本 节 中 将 RGA，SSO, PSO 以 及 RGSO 进行 比较 测试 ， 分 


析 RGSO 的 性 能 。 实 验 中 的 迭代 总 数 设 置 为 


1000 次 ， 每 种 算 


法 分 别 独立 运行 30 次 , 计算 均值 与 方差 。 每 次 测试 当 达 到 最 大 


帮 代 次 数 或 结果 误差 值 小 于 10 取 时 ， 实 验 将 终止 。 


D=60 的 结果 比较 如 表 
示 。 各 算法 在 D=60 上 的 排名 如 图 
如 图 5 所 示 。 


3 所 示 ,D=100 的 结果 比较 如 表 4 所 
4 所 示 ， 在 D=100 上 的 排名 


从 表 3 中 可 以 看 到 ， 在 D=60 时 ，RGSO # fl. f3. f5 以 
及 位 下 都 有 十 分 有 效 的 结果 , 在 其 余 基 准 函 数 下 也 有 较 好 的 结 
果 。 从 表 4 中 可 以 看 到 ， 在 D=100 时 ，RGSO Æ fl, B. f5 以 
及 位 下 都 有 十 分 有 效 的 结果 , 在 其 余 基 准 函 数 下 也 有 较 好 的 结 
Ko 从 图 4 和 图 5 中 可 以 看 到 , 相 较 于 其 他 算法 , RGSO E fl, 
2. 8. f5, f6, f7 上 有 着 更 加 精确 的 结果 。 


在 大 部 分 基准 函数 下 ，RGSO 


了 着 良好 的 效果 且 算法 表现 


较 稳 定 。 但 在 某 些 多 模 函 数 下 ，RGSO 的 表现 


不 十 分 出 色 。 


经 分 析 ， 原 因 如 下 : a)RGSO 的 理论 基础 是 由 一 颗 种 子 发 育 为 


整个 种 群 ， 因 此 ， 初 始 种 子 位 置 的 选取 一 定 程度 上 影响 了 最 终 
算法 的 和 迭代 精度 与 速度 ;pb) 分 支 数 的 选取 对 结果 也 有 影响 。 
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表 2 基准 测试 函数 表 
基准 函数 函数 方程 维度 D) 取 值 范围 ”最 小 值 
D 
Sphere HORDI 60&100 [-100, 100] 0 
isl 
D-1 3 , 
Rosenbrock fr) = J 00004? =x)? + —,)") 60&100 — [-30, 30] 0 
isl 
D D 
Schwefel 2.22 Rœ =F s+ [ [|x 60&100 [-10,10] 0 
i=l isl 
D 
Rastrigin Fœ = > a —10cos(2zx,) + 10) 60&100 [-10,10] 0 
isl 
£ 2 
SumSquares AO= i 60&100 [-10, 10] 0 
i=l 
fi or 
F(x) = 20 + exp- 20exp(-0.2 aoe ) 
Ackley 7 = 60&100 [-32, 32] 0 
~exp(+cos(2zx,)) 
DI 
: _ 1 Aam PEA S 
Griewank A@)= F952) qprogp+ 60&100 [-600, 600] 0 
D 
F4 f = 418.9829n + X =x, sinl |x 60&100  [-100, 100] 0 
i=l 
D p D 
Salomon 万 =-cosCz S oLD 37 +1 60&100  [-100, 100] 0 
isl isl 
D 
fio=0.1(sin’ (3zx) + >) (x, -D° (1+ sin’ (Bax, +1) 
Levy a 608&100 [-10, 10] 0 
+(x, — 1) (1+ sin’ (27x, )) 
#23 RGA, SSO, PSO, RGSO 的 结果 比较 (D=60) #4 RGA, SSO, PSO, RGSO 的 结果 比较 (D=100) 
Function RGA SSO PSO RGSO Function RGA SSO PSO RGSO 
Mean 0.0056 1.0565e+003 1.8134e+004 3.2030e-005 Mean 0.5357 68.9464  4.8886e+004 0.0257 
Sphere Sphere 
Std 0.0023 4.4693e+003 1.3196e+004 1.1855e-005 Std 0.1439 13.6176  2.1782e+004 0.0077 
Mean 279.1657 117.1419 8.0333e+006 69.7145 Mean 1.0521e+004 690.8419 1.4017e+008 136.0489 
Rosenbrock Rosenbrock 
Std 111.7336 41.1482  2.4626e+007 28.6350 Std 2.7700e+004 142.5517  1.0982e+008 40.9249 
Mean 132.1208 374.4894 604.1314 0.0579 Mean 718.5855 689.6575 861.8725 6.289295 
Schwefel 2.22 Schwefel 2.22 
Std 124.1313 211.5519 372.3285 0.0020 Std 164.5298 210.7466 109.4321 2.183518 
Mean 281.6219 291.8714 363.0958 303.9969 Mean 595.8592 645.8945 760.3829 550.0200 
Rastrigin Rastrigin 
Std 43.0766 55.9090 62.8330 27.0229 Std 83.5125 70.6721 57.1647 51.4200 
Mean 0.2279 0.9511 8.5148  4.5992e-005 Mean 12.1975 37.4532 182.6591 0.0672 
SumSquares SumSquares 
Std 0.1898 0.8232 3.5358  4.0291e-005 Std 4.7167 15.5354 83.5327 0.0209 
Mean 4.2865 4.8903 3.7933 2.6506 Mean 5.5952 6.8521 5.8084 3.3317 
Ackley Ackley 
Std 0.4581 0.8484 0.4963 0.5124 Std 0.6533 0.7191 0.6568 0.5793 
Mean 8.7372e-004 0.3704 1.0331  5.8924e-005 Mean 0.0685 1.0633 1.6492 0.0211 
Griewank Griewank 
Std 1.0842e-004 0.0664 0.0306 2.5941e-005 Std 0.0166 0.0386 0.1578 0.0077 
Mean 10.5784 50.3823 103.0523 1.0989e-003 Mean 862.0364 6.1934e+003 4.5236e+006 1.0998 
F4 F4 
Std 6.7566 47.8534 82.0002 0.5067e-003 Std 696.8623 1.1786e+003 2.0902e+006 0.7870 
Mean 941.5325 1.6509e+003 3.6255e+003 1.9854 Mean 6.1263e+004 3.8265e+006 6.5956e+008 55.5455 
Salomon Salomon 
Std 865.2547 1.4654e+003 2.1242e+003 0.8965 Std 4.9462e+004 1.9023e+006 3.7548e+008 29.8974 
Mean 8.4367e-003 0.7243 10.5673  3.8652e-004 Mean 68.6240 993.8613 3.2852e+004 0.2133 
Levy Levy 
Std 5.9672e-003 0.1265 8.1839  2.9524e-004 Std 55.6724 530.0251 2.8090e+004 0.0954 
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5 
mRGA #SSO mPSO mRGSO 
4 
x 
$ 3 
ik 
ak 
xa 2 
1 
0 
fl f2 f3 f4 f5 f6 f7 
D=60 
图 4 各 算法 在 D=60 上 的 排名 
5 
ERGA mSSO npPSO mRGSO 
4 
x 
x 
Ia 
xt 
1 
0 
fl f2 f3 f4 f5 f6 f7 
D=100 
图 5 各 算法 在 D=100 上 的 排名 
5 ”结束 语 


受 启发 于 植物 根系 的 生长 过 程 ， 本 文 在 已 有 根 生 算法 的 基 
础 上 作出 改进 , 提出 了 根 生 群 优化 算法 (RGSO) ,用 于 解决 优 
化 问题 。 通 过 构建 根系 群 模 型 ， 描 述 了 根系 生长 过 程 中 根系 总 
体 的 E 长 模型 和 引入 
SVDD 构建 土壤 养分 分 布 模型 阐述 了 不 同 分 类 的 根 尖 的 生长 方 
式 ,使 得 种 群 内 部 的 类 别 划分 更 加 合理 , 也 加 快 了 算法 的 收敛 速 
度 和 精确 性 。 从 理论 上 证 明了 算法 的 全 局 收敛 性 。 在 实验 中 调 
节 参 数 ， 确 定 了 使 算法 达到 最 佳 效 果 的 参数 值 。 与 当前 最 先进 
的 其 他 算法 的 对 比 ， 如 PSO，SSO，RGA 等 ， 结 果 显 示 RGSO 
具有 较 高 的 有 效 性 以 及 全 局 收敛 性 。 
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